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Abstract. Let A be a self-adjoint operator on a Hilbert space Sj. Assume that the spectrum of 
A consists of two disjoint components cTq and CTi . Let V be a bounded operator on j^, off-diagonal 
and y-self-adjoint with respect to the orthogonal decomposition = ioo ffi^l where S)o and S)i 
are the spectral subspaces of A associated with the spectral sets C7o and CTi , respectively. We find 
(optimal) conditions on V guaranteeing that the perturbed operator L = A + V is similar to a self- 
adjoint operator. Moreover, we prove a number of (sharp) norm bounds on the variation of the 
spectral subspaces of A under the perturbation V. Some of the results obtained are reformulated 
in terms of the Krein space theory. As an example, the quantum harmonic oscillator under a 
J^^'-symmetric perturbation is discussed. 



1. INTRODUCTION 

Let A be a (possibly unbounded) self-adjoint operator on a Hilbert space Sj. Assume that V 
is a bounded operator on ^. It is well known that in such a case the spectrum of the perturbed 
operator L = A + V lies in the closed ||y ||-neighborhood of the spectrum of A even if V is non- 
self-adjoint. Thus, if the spectrum of A consists of two disjoint components Gq and Oi, that 
is, if 

spec(A) = aoUai and dist(ao,ai) = <i > 0, (1.1) 

then the perturbation V with a sufficiently small norm does not close the gaps between Gq and 
d in C. This allows one to think of the corresponding disjoint spectral components Oq and a[ 
of the perturbed operator L = A + V as a result of the perturbation of the spectral sets Oo and Oi , 
respectively. 

Assuming dl.ll ). by E^(ao) and E^(ai) we denote the spectral projections of A associated 
with the disjoint Borel sets Gq and Oi, and by and Sj\ the respective spectral subspaces, 
^0 = RanEA(cJo) and = RanE^(ai). If there is a possibility to associate with the disjoint 
spectral sets Gq and a[ the corresponding spectral subspaces of the perturbed (non-self-adjoint) 
operator L = A + V, we denote them by Sj'q and Sj[ . In particular, if one of the sets cTq and (j[ is 
bounded, this can easily be done by using the Riesz projections (see, e.g. |24, Sec. III.4]). 

In the present note we aie mainly concerned with bounded perturbations V that possess the 
property 

V*=JVJ, (1.2) 
where 7 is a self-adjoint involution on Sj given by 

J=EA{oo)-EA{ai). (1.3) 
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Operators V with the property (11.21 ) are called 7-self-adjoint. 

A bounded perturbation V is called diagonal with respect to the orthogonal decomposition 
^ = ^0 © ^1 if it commutes with the involution J, VJ = JV. If V anticommutes with J, i.e. VJ = 
—JV, then V is said to be off-diagonal. Clearly, any bounded V can be represented as the sum 
y = Vdiag + Voff of the diagonal, Vdiag, and off-diagonal, Voff, terms. The spectral subspaces 
and ^1 remain invariant under A + Vdiag while adding a non-zero Voff does break the invariance 
of ^0 and Thus, the core of the perturbation theory for spectral subspaces is in the study 
of their variation under off-diagonal perturbations (cf. [25]). This is the reason why we add to 
the hypothesis (11.21 ) another basic assumption, namely that all the perturbations V involved are 
off-diagonal with respect to the decomposition ^ = © ^ i ■ 

We recall that if an off-diagonal perturbation V is self-adjoint in the usual sense, that is, 
V* = V, then the condition 

n n d 

ll^ll<2 (1-4) 

ensuring the existence of gaps between the perturbed spectral sets and a[ may be essentially 
relaxed. Generically, if no assumptions on the mutual position of the initial spectral sets Oq and 
Gi are made except (II. lb . the sets Oq and o[ remain disjoint for any off-diagonal self-adjoint V 

satisfying the bound ||V|| < (see ||27l Theorem 1 (i)]). If, in addition to (II. lb . it is known 
that one of the sets (Jq and Oi lies in a finite gap of the other set then this bound may be relaxed 
further: for the perturbed sets Oq and o[ to be disjoint it only suffices to require that ||V|| < V2d 
(see llTl Theorem 2 (i)]; cf. ||26l Remark 3.3]). Finally, if the sets do and ai are subordinated, 
say supao < infai, then no requirements on ||V|| are needed at all: the interval (supao,inf ai) 
belongs to the resolvent set of the perturbed operator L = A + V for any bounded off-diagonal 
self-adjoint V (see ||2l[T7l|33l; cf. ll28l ) and even for some off-diagonal unbounded symmetric V 
(see BTI Theorem 1]). It is easily seen from Example |5 . 5 I below that in the case of 7-self-adjoint 
off-diagonal perturbations the condition (11.41 ) ensuring the disjointness of the perturbed spectral 
sets Gq and g[ can be relaxed for none of the above dispositions of the initial spectral sets Gq 
and G\. 

Assuming that V is a bounded 7-self-adjoint off-diagonal perturbation of the (possibly unbo- 
unded) self-adjoint operator A we address the following questions: 

(i) Does the spectrum of the perturbed operator L = A + V remain real under conditions 
(fTTT) and (O)? 

(ii) If yes, is it then true that L is similar to a self-adjoint operator? 

(iii) What are the (sharp) bounds on variation of the spectral subspaces associated with the 
spectral sets Gq and Gi as well as on the variation of these sets themselves? 

In our answers to the above questions we distinguish two cases: 

(G) the generic case where no assumptions on the mutual positions of the spectral sets Gq 

and Gi are done except for the disjointness assumption dl.ll ): 
(S) the particular case where the sets Gq and Gi are either subordinated, e.g. sup Gq < inf Gi, 

or one of these sets lies in a finite gap of the other set, say Gq lies in a finite gap of Gi . 

We have to underline that this distinction is quite different from the one that arises when the 
perturbations V are self-adjoint in the usual sense: the case (S) now combines the two spectral 
dispositions that should be treated separately if V were self-adjoint (see r7l [T7ll27ll4TI '). 

Our answers to the questions (i) and (ii) are complete and positive in the case (S). In this case 
the spectrum of the perturbed operator L = A + V does remain real for any off-diagonal 7-self- 
adjoint V satisfying the bound ||V|| <d/2. Moreover, the operator L turns out to be similar to a 
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self-adjoint operator whenever the strict inequality (11.41 ) holds. These results combined in The- 
orem [53] (ii) below (see also Remark l5.13l ) represent an extension of similar results previously 
known due to [1] and [37] for the spectral dispositions with subordinated Oq and Oi. 

By using the results of [32^, '49^, we give a positive answer to the question (i) also in the 
generic case (G) provided that the unperturbed operator A is bounded (see Theorem 15. 121 ). For 
A unbounded, we prove that in case (G) the spectrum of L = A + V for sure is purely real if V 
satisfies a stronger bound ||y || <d/7i. The strict bound \\V\\ <d/n guarantees, in addition, that 
L is similar to a self-adjoint operator (see Theorem 15. 8l (i)). The question whether this is true 
for d/n < \\V\\ <d/2 remains an open problem. 

We answer the question (iii) by using the concept of the operator angle between two sub- 
spaces (for discussion of this notion and references see, e.g., [25]). Recall that if 9Jt and are 
subspaces of a Hilbert space, the operator angle 0(9Jt,OT) between Tl and measured relative 
to the subspace dJl is introduced by the following formula |[26l : 



e{m,m)=arcsm^Jln-PwP<n\^, (1-5) 

where Itjji denotes the identity operator on DJl and Ptjji and P<xi stand for the orthogonal projections 
onto Tl and respectively. 
Set 

2d/n, case (G), ^ 
d, case(S), 

and assume that ||y || < 5 /2. Since in both the cases (G) and (S) under this assumption we have 
got the positive answer to the question (ii), one can easily identify the spectral subspaces S)q and 
Sj\ of L associated with the corresponding perturbed spectral sets Gq and a[ (cf. Lemma [5^ . 
Let &j = @{S^j,9)'j), J = 0, 1, be the operator angle between the unperturbed spectral subspace 
S)j and the perturbed one, S^'j. Our main result (presented in Theorem |5.8l ) regarding the operator 
angles ©o and ©i is that under condition ||V|| < 8/2 the following bound holds: 

fl 2||y||\ 

tan 0^- < tanh ( - arctanh j , j = 0, 1 , (1 .7) 

which means, in particular, that 0y < |, j = 0, 1. Theorem [5^ also gives the bounds on location 
of the perturbed spectral sets and o[ (see formulas (15.191 )) . 

In the case (S) the bounds on a,, and a[ as well as the bounds (11.71 ) are optimal (see Re- 
mark l5.10l) . Inequalities (11.71 ) resemble the sharp norm estimate for the operator angle between 
perturbed and unperturbed spectral subspaces from the celebrated Davis-Kahan tan 20 Theo- 
rem (see fTTl . p. 11; cf. |[28l Theorem 2.4] and BT[ Theorem 1]). Recall that the latter theorem 
serves for the case where the unperturbed spectral subsets Oq and 0\ are subordinated and the 
off-diagonal perturbation V is self-adjoint. The difference is that the usual tangent of the Davis- 
Kahan tan 20 Theorem is replaced on the right-hand side of (11.71 ) by the hyperbolic one. Another 
distinction is that the bound ( 11.71 ) holds not only for the subordinated spectral sets Oq and ai but 
also for the disposition where one of these sets lies in a finite gap of the other set and thus Oq 
and 0\ are not subordinated. 

The results obtained are of particular interest for the theory of operators on Krein spaces |[9l . 
The reason for this is that introducing an indefinite inner product [x,y\ = {Jx,y), x,y G S), instead 
of the initial inner product (•,•), turns 9) into a Krein space. The operators V and L = A + V being 
7-self-adjoint on f) appear to be self-adjoint operators on the newly introduced Krein space ^. 
Under the condition ||y || < 5/2 in both cases (G) and (S) we establish that the perturbed spectral 
subspaces and are mutually orthogonal with respect to the inner product [•,•]. Moreover, 
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these subspaces are maximal uniformly positive and maximal uniformly negative, respectively 
(see Remark lS.llI ). The restrictions of L onto S^'q and are .^-unitary equivalent to self-adjoint 
operators on and ^i, respectively. This extends similar results previously known from HI 
and ||371 for the case where the spectral sets Oq and 0\ are subordinated. 

Another motivation for the present paper is in the spectral analysis of non-self-adjoint Schro- 
dinger operators that involve the so-called ^-symmetric potentials. Starting from the pio- 
neering works flT, T2], these potentials attracted considerable attention because of their prop- 
erty to produce, in some cases, purely real spectra (see, e.g., |3, 4(Ll0,,29ii39,i50J). The local 
^=^-symmetric potentials appear to be /-self-adjoint with respect to the space parity operator 
^ (see, e.g., ll32l[391 ). allowing for an embedding the problem into the context of the spectral 
theory for 7-self-adjoint perturbations (this also means that the ,^^=;7-symmetric perturbations 
may be studied within the framework of the Krein space theory Bl [32ll48 l). 

The main tool we use in our analysis is a reduction of the problems (i)-(iii) to the study of 
the operator Riccati equation 

KAq +AiK + KBK = -B* 

associated with the representation of the perturbed operator L = A + V in the 2x2 block matrix 
form 

where Aq = A|^^, Ai = A|^ , and B = V |^ . Assuming (11.61 ). we prove that the Riccati equation 
has a bounded solution K for any B such that ||B|| < 5/2. The key statement is that the perturbed 
spectral subspaces and are the graphs of the operators K and K* , respectively, which then 
allows us to derive the bounds (11.71) . 

The plan of the paper is as follows. In Section |2] we give necessary definitions and present 
some basic results on the operator Riccati equations associated with a class of unbounded non- 
self-adjoint 2x2 block operator matrices. Section[3]is devoted to the related Sylvester equations. 
In Section |4] we prove a number of existence and uniqueness results for the operator Riccati 
equations. In Section [5] we consider 7-self-adjoint perturbations and find conditions on their 
norm guaranteeing the reality of the resulting spectrum. In this section we also prove the bound 
(11.71) on the variation of the spectral subspaces and discuss the embedding of the problem into the 
context of the Krein space theory. Finally, in Section [6] we apply some of the results obtained to 
a quantum-mechanical Hamiltonian describing the harmonic oscillator under a ^^^-symmetric 
perturbation. 

We conclude the introduction with the description of some more notations that are used thro- 
ughout the paper. By a subspace we always understand a closed linear subset of a Hilbert space. 
The identity operator on a subspace (or on the whole Hilbert space) 9Jt is denoted by /jjf. If no 
confusion arises, the index 971 may be omitted in this notation. The Banach space of bounded 
linear operators from a Hilbert space OJt to a Hilbert space is denoted by =^(9Jt, OT). For 
we use a shortened notation =^(9Jt). By 9Jt©9T we will understand the orthogonal 
sum of two Hilbert spaces (or orthogonal subspaces) SOT and OT. By (^^(SDT, DT), < r < oo, we 
denote the closed ball in =^(971,91), having radius r and being centered at zero, that is, 

^,(971, = {Ke ^{m,m) I \\K\\ < r}. 

If it so happens that r = +oo, by ^c»(9Jt,DT) we will understand the whole space ^(9Jt,9T). 
The notation conv(a) is used for the convex hull of a Borel set a C M. By Or{^), r > 0, we 
denote the closed r-neighborhood of a Borel set Q. in the complex plane C, i.e. Or{^) = {z G 
C| dist(z,a) < r}. 
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2. Operator Riccati equation 

We start by recalling the concepts of weak, strong, and operator solutions to the operator 
Riccati equation (see ||5l|6l). 

Definition 2.1. Assume that Aq and Ai are possibly unbounded densely defined closed operators 
on the Hilbert spaces and i^i , respectively. Let B and C be bounded operators from to 
and from to i^i, respectively. 

A bounded operator K € =^(^0,^1) is said to be a weak solution of the Riccati equation 

KAq-AiK + KBK = C (2.1) 

if 

{KAQX,y) - [Kx,A\y) + {KBKx,y) = {Cx,y) 
for all X £ Dom(Ao) and y £ Dom(Aj). 
A bounded operator K G ^(ijo, ^1 ) is called a strong solution of the Riccati equation ( 12.11 ) if 

Ran(^|Dom(Ao)) CDom(Ai) (2.2) 

and 

KAqx-AiKx + KBKx = Cx for all ;ceDom(Ao). (2.3) 
Finally, K G ^(i5o,^i) is said to be an operator solution of the Riccati equation (12.11 ) if 

Ran(^) CDom(Ai), (2.4) 
the operator KAq is bounded on Dom{KAo) = Dom(A()), and the equality 



KAq-AiK + KBK = C (2.5) 



holds as an operator equality, where KAq denotes the closure of KAq. 
Remark 2.2. We will call the equation 

XA*i -A*oX-XB*X = -C* (2.6) 

the adjoint of the operator Riccati equation (12.11 ). It immediately follows from the definition 
that an operator K G ^{S)o,Sj\) is a weak solution to the Riccati equation (12.11 ) if and only if 
the adjoint of K, X = K* , is a. weak solution to the adjoint equation ( 12.61 ). 

Clearly, any operator solution K G ^(^0,^1) to the Riccati equation (12.11 ) is automatically a 
strong solution. Similarly, any strong solution is also a weak solution. But, in fact, by a result of 
one does not need to distinguish between weak and strong solutions to the Riccati equation 
(12.11 ). This is seen from the following statement. 



Lemma 2.3 ([6|, Lemma 5.2). Let Aq and Ai be densely defined possibly unbounded closed 
operators on the Hilbert spaces ?)() and ?)\, respectively, and B £ C £ ^{S^q,?)i). If 

K £ ^{:^Q,S)i) is a weak solution of the Riccati equation (12.11 ) then K is also a strong solution 



If the operators Aq, Ai, B, and C are as in Definition 12.11 then a 2 x 2 operator block matrix 
Ao B 



., Dom(L) =Dom(Ao)eDom(Ai), (2.7) 

is a densely defined and possibly unbounded closed operator on the Hilbert space 

^ = ^00^1. (2.8) 
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The operator L will often be viewed as the result of the perturbation of the block diagonal matrix 
A = diag(Ao,Ai), Dom(A) = Dom(Ao) ©Dom(Ai), (2.9) 
by the off-diagonal bounded perturbation 

v = {ll)- a.io) 

The operator Riccati equation (12.11 ) and the block operator matrix L are usually said to be 
associated to each other. Surely, one can also associate with the matrix L another operator 
Riccati equation, 

K'A 1 - AqK' + K'CK' = B, (2.11) 

assuming that a solution K' (if it exists) should be a bounded operator from to ^o- 

It is well known that the solutions to the Riccati equations (12.11) and (12.111) determine invariant 
subspaces for the operator matrix L (see, e.g., |0 for the case where the matrix L is self-adjoint 
or on for the case of a non-self-adjoint L). These subspaces have the form of the graphs 

'i^{K) = {x e SjQ® S)i\x = X(jBKx(j for some xq G ^o} (2.12) 

and 

^{K') = {xeS3o®S)i\x = K'xi®xi for some ;ci G (2.13) 

of the corresponding (bounded) solutions K and K'. Notice that the subspaces of the form 
(12.121) and (12.131) are usually called the graph subspaces associated with the operators K and K', 
respectively, while K and K' themselves are called the angular operators. Usage of the latter 
term is explained, in particular, by the fact that if a subspace (S C ^ is a graph 6 = 'i^{K) of a 
bounded linear operator K from a subspace dJl to its orthogonal complement 9Jl^, 9Jl^ = ^ 9Jt, 
then the following equality holds (see |25|; cf. ||T3 and EB): 

|/:| = tan 0(971,0), (2.14) 

where |^| is the absolute value of K, \K\ = \/K*K, and 0(S!Jt, 6) the operator angle (11.51 ) be- 
tween the subspaces 9Jt and 0. 

The precise statement relating solutions of the Riccati equations (12.11 ) and (12.111 ) to invariant 
subspaces of the operator matrix (12.71 ) is as follows. 

Lemma 2.4. Let the entries Aq, A\, B, and C be as in Definition \2.1\ and let a 2x2 block 
operator matrix L be given by (12.71 ). Then the graph ^{K) of a bounded operator K from 
S^Q to ^1 satisfying ( 12.21 ) is an invariant subspace for the operator matrix L if and only if K 
is a strong solution to the operator Riccati equation (12.11 ). Similarly, the graph ^{K') of an 
operator K' £ .^{S)i,?)o) such that Ran(/r'|£)oni(^|)) C Dom(Ao) is an invariant subspace for L 
if and only if this operator is a strong solution to the Riccati equation (12.111 ). 

The proof of this lemma is straightforward and follows the same line as the proof of the 
corresponding part in [T, Lemma 5.3]. Thus, we omit it. 

The next assertion contains two useful identities involving the strong solutions to the Riccati 
equations (12.11 ) and (12.111 ). 

Lemma 2.5. Let the entries Aq, A\, B, and C be as in Definition \2.1\ Assume that operators K G 
^{f)o,Sji) and K' G ^{S^i,Sjq) are strong solutions to equations (12.11 ) and (12.111) . respectively. 
Then 

Ran(^'^|Dom(Ao)) CDom(Ao), Ran(/:/:'|Do^(^,)) C Dom(Ai), (2.15) 
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and 

{I - K'K) (Ao + BK)x = (Ao - K'C) (/ - K'K)x for all x G Dom(Ao) , (2. 16) 
{I-KK'){Ai+CK')y = {Ai-KB){I-KK')y /or a/Z 3; e Dom(Ai). (2.17) 

Proof. The inclusions (12.151 ) follow immediately from the definition of a strong solution to the 
operator Riccati equation (see condition (12.21 )). 

Let X G Dom(Ao). Taking into account the first of the inclusions (12.151) as well as the inclu- 
sions Ran(A'|Do^(^u)) C Dom(Ai) and Ran(A''|QQ^(^j)) C Dom(Ao) one can write 

(Ao - K'C) (/ - K'K)x = (Ao - K'C)x - {AqK' - K'CK')Kx 

= {A(i-K'C)x-{K'Ai -B)Kx, (2.18) 

by making use of the Riccati equation (12.111 ) itself at the second step. Similarly, 

(/ - K'K) (Ao + BK)x = (Ao + BK)x - K'{KAq + KBK)x 

= {Aq + BK)x-K'{C+AiK)x, (2.19) 

due to the Riccati equation (12.11 ). Comparing (12.181) and (12.191) we arrive at the identity (12.161 ). 
Identity (12.171 ) is proven analogously. □ 

We will also need the following auxiliary lemma. 

Lemma 2.6. Suppose that the operators K G ^{S)q,Sji) and K' E ,'3§(Sji,Sj()) are such that the 
2x2 operator block matrix 

^={k f) ^2.20) 

considered on = S)o(B S)i is boundedly invertible, i.e. the inverse operator W^^ exists and is 
bounded. Then the graphs '^{K) and ^S{K') of the operators K and K' are linearly independent 
subspaces ofSj and 

S^ = ^{K)\^{K'), (2.21) 
where the sign "+" denotes the direct sum of two subspaces. 

Proof. The existence and boundedness of imply that equation Wx = yis uniquely solvable 
for any y ^ S). This means that there are unique xq G and unique xi ^ Sji such that y = 
xqBKxo + K'xi 0X1 and hence ^ C '^{K) +'^{K'). Since both '^{K) and ?^(/^') are subspaces 
of S), the inclusion turns into equality, S) = "^{K) +'^{K'). The linear independence of '^{K) 
and W{K') follows from the fact that equation Wx = has only the trivial solution x = 0. □ 

Remark 2.7. It is well known that the following three statements are equivalent. 

(i) The operator matrix (12.201 ) is boundedly invertible. 

(ii) The inverse (/ — KK')^^ exists and is bounded. 

(iii) The inverse {I — K'K)^^ exists and is bounded 

For a proof of this assertion see, e.g. 1221 Theorem 1.1. and Lemma 2.1] where even a Banach- 
space case of 2 x 2 block operator matrices of the form (12.201) with unbounded entries K and K' 
has been studied. 

Remark 2.8. The inverse of the operator W is explicitly written as 

{I-K'K)~' -K'iI-KK')-^\ 
~ \-K{I-K'K)-^ {I-KK')-^ J 
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{I-K'K)-' -{I-K'K)-'K'\ 

-{1-KK')-^K {I-KK')-\ ) ^ ' 

The (oblique) projections Qi^i^k) ^^id Q'S(k') oi^to the graph subspaces 'SiK) and "^{K') along 
the corresponding complementary graph subspaces '^{K') and ^(K) are given by 

= (j^) {i-K'K)-' (7 -K') and Q^^^^') = (f ) {I - KK'y' {-K l) , (2.23) 
respectively. 

Corollary 2.9. Assume the hypothesis of Lemma \2.4\ Suppose that K G =^(^0)^i) '^^d K' G 
•^(^1,^0) fl'''^ strong solutions to the Riccati equations ( 12.11 ) and ( 12.1 11 ). respectively. Assume, 
in addition, that the 2x2 operator block matrix W formed with these solutions according to 
(12.201) is a boundedly invertible operator on S) = iDo ©^i- Then: 

(i) The operator L is similar to a block diagonal operator matrix Z = diag(Zo,Zi), 

L = WZW-\ (2.24) 
where Zq and Z\ are operators on 9)q and Sji, respectively, given by 

Zo=Ao + BK, Dom(Zo) =Dom(Ao), (2.25) 

Zi=Ai+CK', Dom(Zi) =Dom(Ai). (2.26) 

(ii) The Hilbert space Sj splits into the direct sum S) = Sj\ of the graph subspaces 
Sjq = ^{K) and 9)\ = ^{K') that are invariant under L. The restrictions L\^i^ andL\^i^ 
ofL onto S^'q and ^\ are similar to the operators Zq and Z\, 

Wo-1l|j5/Wo=Zo and IVf ^Llj^.m = Zi, (2.27) 

where the entries Wb : ^0 ^ ^0 Wi : ^ 1 ^ ^ j correspond to the respective columns 
of the block operator matrix W , 

Wb^co = -^0, -^0 G ^0, and Y^\X\ = ^ xi , xi G ^1. (2.28) 

Proof. First, one verifies by inspection that LW = WZ taking to account that K and K' are the 
strong solutions to the Riccati equations (12.11 ) and (12.111 ). respectively. The remaining statements 
immediately follow from Lemma 124] combined with Lemma l26l □ 

Remark 2.10. The similarity (12.241) of the operators L and Z implies that the spectrum of L 
coincides with the union of the spectra of Zq and Z\, that is, spec(L) = spec(Zo) U spec(Zi). 

3. Operator Sylvester equation 

Along with the Riccati equation (12.1b we need to consider the operator Sylvester equation 

XAq-AxX = Y (3.1) 

assuming that the entries Aq and A\ are as in Definition 12. II and Y G =^(i5o)^i)- The Sylvester 
equation is a particular (linear) case of the Riccati equation and its weak, strong, and operator 
solutions X G =^(^0)^i) are understood in the same way as in the above definition. Further- 
more, by Lemma 123] (cf. |8, Lemma 1.3]) one does not need to distinguish between the weak 
and strong solutions to (13.11 ). 

Because of its importance for various areas of mathematics there is an enormous Uterature on 
the Sylvester equation (for a review and many references see paper lfT4l '). With equation (13.11 ) 
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one often associates the Sylvester operator S defined on the Banach space =^(^0,^1) by the 
left-hand side of (1311) : 

S{X)=XAo-AiX (3.2) 

with domain 

Dom(S) = {x E .^{^0,^1) I Ran(X|p^^(^^p C Dom(Ai)} . (3.3) 

Clearly, the Sylvester equation (13.11) has a unique solution X e Dom(S) if and only if 
spec(S). It is known that in general the spectrum of S is larger than the (numerical) difference 
between the spectra of Aq and Ai. More precisely, provided that spec(Ao) 7^ C or spec (A 1) / C 
always the following inclusion holds HI : 

spec(Ao) — spec(Ai) C spec(S), (3.4) 

where we use the notation £ — A = {z — | z G Z, G A} for the numerical difference between 
two Borel subsets £ and A of the complex plane C. The opposite inclusion in (13.41 ) may fail to 
hold if both operators Aq and Ai are unbounded. The corresponding example was first given 
by V. Q. Phong [43] for the Sylvester equation (13.11 ) where one of the entries Aq and A 1 is an 
operator on a Banach (but not Hilbert) space. An example where both Aq and Ai are operators on 
Hilbert spaces and spec(S) ^ spec(Ao) — spec(Ai) may be found in IH Example 6.2]. Equality 

spec(S) = spec(Ao) — spec(Ai) (3.5) 

holds if both Aq and Ai are bounded operators. This result is due to G. Lumer and M. Rosenblum 
|[34l . Equality (13.51) also holds if only one of the entries Aq and Ai is a bounded operator [8]. In 
this case (13.51 ) implies that if the spectra Aq and Ai are disjoint then spec(S) and hence the 
operator S is boundedly invertible. Moreover, a unique solution of the Sylvester equation (13.11 ) 
admits an "explicit" representation in the form a contour integral. 

Lemma 3.1. Let Aq be a possibly unbounded densely defined closed operator on the Hilbert 
space 9)0 and A{ a bounded operator on the Hilbert space f)\ such that 

spec(Ao) n spec (A 1) = 

and Y G =^(^0)^i)- Then the Sylvester equation (13.11) has a unique operator solution 

X = ^ fdziAi-z)-'Y{Ao-z)-\ (3.6) 

where 7 is a union of closed contours in C with total winding numbers around spec(Ao) and 1 
around spec (A 1 ) and the integral converges in the norm operator topology. 

Corollary 3.2. Under the hypothesis of Lemma \3.1\ the norm of the inverse of the Sylvester 
operator S may be estimated as 

||S-i||<(27r)-i|7|sup||(Ao-z)-^||||(Ai-z)-i||, 

zey 

where \ y\ denotes the length of the contour 7 in (13.61) . 

The result of Lemma 13.11 may be attributed to M. G. Krein who lectured on the operator 
Sylvester equation in late 1940s (see 11141 '). Later, it was independently obtained by M. Rosen- 
blum [45 1 . 

As for the Sylvester operator (13.21 ) with both unbounded entries Aq and Ai, we have an im- 
portant result which is due to W. Arendt, F. Rabiger, and A. Sourour (see (HI Theorem 4.1 and 
Corollary 5.4]). 
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Theorem 3.3 (HI). Let Aq andA\ be closed densely defined operators on the HUbert spaces 
and ^1, respectively. Assume that one (or both) of the following holds (hold) true: 

(i) A() and (— Ai) are generators of eventually norm continuous C{)- semigroups; 

(ii) Aq and (— Ai) are generators of Co-semigroups one of which is holomorphic. 
Then the spectrum of the Sylvester operator (13.21 ) is given by (13.51 ). 

Recall that an operator H on the Hilbert space 9Jt is said to be m-dissipative if it is closed 
and both the spectrum and numerical range of H are contained in the left half -plane {z G C | 
Imz < 0}. The Lumer-Phillips theorem asserts (see, e.g., |[T9l Section II.3.b]; cf. |[20l Theo- 
rem B.21]) that a Co-semigroup on OK is a contraction semigroup if and only if its generator is 
an m-dissipative operator. The next statement represents a generalization of a well known result 
by E. Heinz ( [23. Satz 5]) to the case of unbounded operators. Notice that the exponential e^', 
f > 0, is understood below as the corresponding element of the strongly continuous contraction 
semigroup generated by an (unbounded) m-dissipative operator H. 

Theorem 3.4. Let Aq + |/ and — Ai + f /, 5 > 0, be m-dissipative operators on the Hilbert 
spaces S^o and Sji, respectively, and Y £ ^{Sjo,f)i). Then the Sylvester equation (13.11 ) has a 
unique weak (and hence unique strong) solution given by 

X = - dte-^"Ye'^°', (3.7) 

JO 

where the integral is understood in the weak operator topology. Moreover, the norm of the 
solution (13.71 ) satisfies the estimate 

ll^ll < tII^'II- (3-8) 
o 

Proof. Under the hypothesis the operators Aq and (— Ai ) are themselves m-dissipative. Let Uo{t) 
and Ui{t), f > 0, be contraction Co-semigroups generated respectively by Ao and (— Ai), that is, 
Uo{t) = e'^o' and Ui{t) = q-"^''. Clearly, 

\\Uo{t)\\ < e"^' and \\Ui{t)\\ < t > 0. (3.9) 

The same bound also holds for the adjoint semigroup Ui{t)*, t > 0, whose generator is the 
m-dissipative operator {—A\). Pick up arbitrary x G ^o and y € Sji and introduce the orbit maps 
t = ^^^0(0-^ ai^d '■ f ^ Cv(0 = Ui{t)*y. By the definition of a strongly continuous 

semigroup, these maps are continuous functions of f G [0,°°). Taking into account the bounds 
(13.91 ) one then concludes that the improper integral 

poo poo 

/ dt{Ui{t)YUoit)x,y)= dt{YUt),Ut)) (3.10) 
Jo JO 

converges and its absolute value is bounded by ||y || Thus, the weak integral on the 

right-hand side of ( 13.71 ) exists and the bound ( 13.81 ) holds. 

Now assume that x G Dom(Ao) and y G Dom(Aj). In this case the orbit maps ^x{t) and l^y{t) 
are continuously differentiable in t and jj^xit) = Uo{t)Aox, j^^y{t) = -Ui{t)*A\y. ForX given 
by (I3.7I ). an elementary computation shows that 

poo 

{XAox,y) - {Xx,Aly) = - ^ dt ((F|^,(0, C>.(0) + (YUt), ^yi^))) 

poo 

= - dti{YUtUy{t)) = {Ym,mi 

J 
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taking into account (13.91 ) in the last step. Since ^x{^) = x and Cv(0) = by Definition 12.11 
this implies that the integral (13.71) is a weak (and hence strong) solution to the Sylvester equa- 
tion (|3l1 ). 

To prove the uniqueness of the weak solution (13.71 ) it is sufficient to show that the homoge- 
neous Sylvester equation XAq —AiX = has the only weak solution X = 0. For a weak solution 
X to this equation we have 

(XA()M,v) - {Xu,A*iv) = for all u £ Dom(Ao) and v G Dom(At). (3.11) 

Take the vectors u and v of the form u = ^x{t) = &^°'x, v = (^y{t) = e^^'^'y, t >0, where the orbit 
maps ^x{t) and i^y{t) correspond to some x G Dom(Ao) and y G Dom(Aj) and hence are both 
continuously differentiable in ? G [0,0°). Notice that the assumption x G Dom(Ao), y G Dom(Aj) 
also implies (see, e.g. 1.19., Lemma 1.3]) that u G Dom(Ao), v G Dom(Ai), and Aqu = ^^£,x{t), 
AjV = — ^^(f). With such a choice of u and v it follows from (13.1 II) that 

jj{X^xit),l^yit)) = whenever x G Dom(Ao) and y G Dom(Ai). 

Hence the function {X^xit),l^y{t)), f > 0, is a constant. Moreover, it equals zero since it vanishes 
as f — > oo. This yields in particular that 

{Xx,y) = (X(^,,(0), Cy(0)) = for all x G Dom(Ao) and y G Dom(A^). 

The latter implies X = 0, which completes the proof. □ 

Remark 3.5. A statement similar to Theorem 13.41 was previously announced without a proof 
in JSl (see E Lemma 2.6]). 

The second important example where a bound of the (13.81) type exists is given in |[T3l Theo- 
rem 3.2]. This example is as follows. 

Theorem 3.6. Assume that the operators Aq and Ai are densely defined and closed. Assume, 
in addition, that there is X G p{A\) such that ||Ao — A || < r and || (A\ — A)^^ || < ('"+ S)^^ for 
some r > and 5 > 0. Then for any Y G c'^{SjQ,f) i ) the unique strong solution X to the Sylvester 
equation (13.11 ) admits the estimate 8\\X\\ < \\Y\\. 

If the operators Aq and Ai are normal then no reference point A is needed and the result is 
stated in a more universal form (see |[T3l Theorem 3.2]). 

Corollary 3.7. Let both Aq and Ai be normal operators such that spec(Ao) is contained in a 
closed disk of radius r, r>0, while spec(Ai) is disjoint from the open disk (with the same center) 
of radius r + 8, 8 > Q. Then for any Y G ^{S)o,f)i) the Sylvester equation (13.11 ) has a unique 
strong solution X and 8\\X\\ < \\Y\\. 

The above two theorems and the corollary give examples where the bounded inverse of the 
Sylvester operator S exists and for the norm of S^' the estimate 5||S^' || < 1 holds with some 
8 > 0. Moreover, this estimate is universal in the sense that it remains valid for any Aq and Ai 
satisfying the corresponding hypotheses. 

4. Existence results for the Riccati equation 

In this section we return to the operator Riccati equation (12.11 ) to prove some sufficient con- 
ditions for its solvability. In their proof we will rely just on the assumption that an estimate Uke 
(13.81 ) holds for the solution of the corresponding Sylvester equation. 
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Theorem 4,1, Let Aq and A\ be possibly unbounded closed densely defined operators on the 
Hilbert spaces Sjq andSji, respectively. Assume that the Sylvester operator S defined on ^(^0,^1) 
by (13.21 ) and (13.31) is boundedly invertible (that is, ^ spec(S)) and 

\\S-'\\<^ (4.1) 

for some 5 > 0. Assume, in addition, the operators B G =^(^1,^0) andC G ^{9)q,S)\) are such 
that the following bound holds: 

v1MM<|- (4.2) 

Then the Riccati equation (12.11) has a unique strong solution in the ball i^5/{2||fi||)(^i)^o). The 
strong solution K satisfies the estimate 

\\K\\ < 1'^'' (4.3) 

| + Vx-ll^llll<^ll 

Proof. If B = then the assertion, including the estimate ( 14.31 ). follows immediately from the 
hypothesis on the invertibility of S on ^{9jo,^i) taking into account the bound ( 14.11 ). 

Suppose that B 7^ 0. In this case the proof is performed by applying Banach's Fixed Point 
Theorem. First, we notice that the bounded invertibility of S on J^{Sjo,S)i) allows us to rewrite 
the Riccati equation (12.11) in the form 

K = F{K) 

where the mapping F : ^{S)o,Sji) Dom(S) is given by 

F{K) = S-\C-KBK). 

By (14.11 ) we have 

<^(||C|| + ||B||||/:f), Kei^{S3o,Sji) (4.4) 

and 

\\F{Ki)-F{K2)\\ < ^\\B\\ {\\Ki\\ + 11^211) 11^1 -K2\\, ^1,^2 G ^iS^o,^i)- (4.5) 
The bound (14.41) implies that F maps the ball (^^(^Oj^i) into itself whenever 

||fi|k^ + ||C|| < r5. (4.6) 

At the same time, from (14.51 ) it follows that F is a strict contraction of the ball ^r(^ij^o) 
whenever 

2\\B\\r<5. (4.7) 

Solving inequahties (14.61) and (14.71 ) one concludes that if the radius r of the ball ^r(^i,^o) is 
within the bounds 

<r<-^, (4.8) 

f + Vf-llfi|lllc|| 

then F is a strictly contractive mapping of the ball i^r(^ii^o) into itself. Applying Banach's 
Fixed Point Theorem one then infers that equation (12.11 ) has a unique solution within any ball 
^r(^i)^o) whenever the radius r satisfies ( 14.81 ). This means that the fixed point is the same for 
all the radii satisfying (14.81) and hence it belongs to the smallest of the balls. This conclusion 
proves the bound (14.31) and completes the whole proof. □ 

Remark 4,2, In (I42l)-(|43]) one may set 5 = ||S"' ||"^ 
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Remark 4.3. By using the hyperbolic tangent function and its inverse the bound (14.31) (for 
B / 0) can be equivalently written in the hypertrigonometric form 



\\C\\ / 1 2x/||B||||C|| \ 

K\\<\ j-jr tanh - arctanh ^ " ^ . (4.9) 



Notice that under condition (14.21 ) we always have 



/ 1 2V||B||||C|| \ 

tanh -arctanh ^ ' < 1. 

\2 5 ) 

Remark 4.4. Fixed-point based approaches to prove the solvability of the operator Riccati equa- 
tion with bounded entries Aq and A\ have been used in many papers (see, e.g., ifTSl . ll42l . 
B6l . HTl). In the case where at least one of the entries Aq and Ai is an unbounded self-adjoint or 
normal operator, a fixed-point approach has been employed in fSl], |l6l|, |l36l, and BO]. Theorem 
14. II represents an extension of the fixed-point existence results obtained in pTl Theorem 3.5] 
and ll42l Theorem 3.1] for the Riccati equation (12.11) with both bounded Aq and Ai to the case 
where the entries Aq and A\ are not necessarily bounded. 

Theorem 4.5. Assume the hypothesis of Theorem \4. 1\ Then the block operator matrix L defined 
by (12.71 ) is block diagonalizable with respect to the direct sum decomposition S) = ^ {K) + ^ {K') 
where K is the unique strong solution to the Riccati equation (12.11 ) within the operator ball 
^5/(2||B||)(^0)^i) <^f^d K' the unique strong solution to the Riccati equation (12.111) within the 
operator ball ^5/(2||c||)(^i,^o)- 

Proof. By Theorem 14. II for K the estimate (14.31 ) holds. By the same theorem for K' we have 

l^'ll < ||B|| ( ^ + \/^-||B||||C|| 1 . (4.10) 




Then the hypothesis ||B||||C|| < 8/2 also implies that ||i^|| ||/^'|| < 1. Hence by Remark [277] the 
operator W in (12.201 ) is boundedly invertible. Applying Corollary [Z9] completes the proof. □ 

Remark 4.6. If, in addition, both operators Aq and A\ are normal then, for r defined by 

l|fi||||C|| , (l 2VI|B||||C||\ 

" = ^/||R|| ||r|| fanh -^rpt^nh ^ " (4.11) 



| + Vf-||B||||C|| 




the spectrum of the block matrix L lies in the closed r-neighborhood of the spectrum of its main- 
diagonal part A = diag(Ao,Ai). That is, dist(z,spec(Ao) U spec(Ai)) < r whenever z G spec(L). 
This immediately follows from the representation (I2.24I) - (I2.26I ) and the bounds (14.31 ) and (14.101 ) 
(see also Remark l2?T0l) . Notice that if S / and C / then r < y^||B||||C|| and hence r < \\V\\ 
taking into account that ||y || = max(||B||, ||C||). 

From now on we assume that the entries Aq and Ai are self-adjoint operators with disjoint 
spectra and thus adopt the following 

Hypothesis 4.7. Let Aq and Ai be (possibly unbounded) self-adjoint operators on the Hilbert 
spaces ^Q and ^\ with domains Dom(Ai) and Dom(Ai), respectively. Assume that the spectra 
of the operators Aq and A\ are disjoint and let 

<i = dist(spec(Ao),spec(Ai)) (>0). (4.12) 
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Hypothesis 14.71 imposes no restrictions on the mutual position of the spectral sets spec(Ao) 
and spec(Ai) except that they are disjoint and separated from each other by a distance d. Some- 
times, however, we will consider particular spectral dispositions described in 

Hypothesis 4.8. Assume Hypothesis \4.7\ Assume, in addition, that either the spectra o/Aq and 
Ai are subordinated, that is, 

sup spec(Ao) < inf spec(Ai) or inf spec(Ao) > sup spec(Ai), (4.13) 

or one of the sets spec(Ao) and spec (A i) lies in a finite gap of the other set, that is, 

conv(spec(Ao)) nspec(Ai) = or spec(Ao) nconv(spec(Ai)) = 0. (4.14) 

Under Hypotheses 14.71 or 148] the bound on the norm of the inverse of the Sylvester operator 
(13.21 ) may be given in terms of the distance d between spec(Ao) and spec(Ai). The following 
result is well known. 

Theorem 4.9. Assume Hvpothesis \4.7\ Let the Sylvester operator 5 be defined by (13.21) and (13.31) . 

(i) Then the inverse of S exists and is bounded. Moreover, the following estimate holds: 

l|S-'ll<^. (4.15) 

(ii) Assume Hvpothesis \4.8\ Then the following stronger inequality holds: 



\\S-'\\<]j. (4.16) 
d 

Remark 4.10. In the generic case (i), where no assumptions on the mutual position of the sets 
spec(Ao) and spec(Ai) are imposed, the existence of a universal constant c such that ||S^' || < — 

has been proven in lfT3l . The proof of the fact that c = ;r/2 is best possible is due to R. McEachin 
||35]| . For more details see 15] Remark 2.8]. As for the particular spectral disposition ( 14.131) . 
the bound (14.161 ) is an immediate corollary to Theorem 13.41 Since any self-adjoint operator 
is simultaneously a normal operator, in the case of the spectral disposition ( 14.141 ) the bound 
(14.161) follows from Corollary 13.71 Sharpness of the bound (14.161 ) in case (ii) is proven by an 
elementary example where the spaces and are one-dimensional, = ^i = C, and the 
entries Aq = ao and Ai = ai are real numbers such that \ai — £?o| = > 0- 

Under the assumption that both the entries Aq and Ai are self-adjoint operators, below we 
present an existence result for the operator Riccati equation ( 12.11 ), which is written directly in 
terms of the distance between the spectra of the entries Aq and Ai (and norms of the operators B 
and C). The result is an immediate corollary to Theorems 14. 1 1 and 14.91 We only notice that the 
role of the quantity 5 in the bounds like (14.11 ) (see inequalities (14.181) and (14.201 ) below) will be 
played by either ^d from (14.151) or d from (14.161 ). 

Theorem 4.11. Assume Hvpothesis \4. 71 

(i) Then for any B G ^{S)i,S)o) andC £ £i^{S)o,S^i) such that 

d 



V\\B\\\\C\\ < - (4.17) 

the Riccati equation (12.11 ) has a unique strong solution K in the ball ^d/{n\\B\\){^OT^i)- 
This solution satisfies the estimate 



\K\\ < , (4.18) 

d , /d^ 

71 



+ ,/S-||B||||C|| 
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(ii) If the conditions of Hypothesis 14.81 also hold then the Riccati equation (12.11 ) has a 
unique strong solution K in the ball ^d/(2||B||) (^0,^i) whenever B G ,!^{Sji,S)o) and 



C G ^{Sjo,S^i) satisfy the bound 



d 



V\\B\\\\C\\<-. (4.19) 
The solution K satisfies the estimate 

\\K\\ < -JE^^^=. (4.20) 

i+VT-iifiiiiicii 

Remark 4.12. The part (i) is a refinement of Tiieorem 3.6 in 1 5 1 tiiat only claimed the existence 
of a weak (but not strong) solution to the Riccati equation (12.11) within the ball ^d/{7t\\B\\){^i:^o)- 
The result of the part (ii) is new. 

Remark 4.13. Let r be given by formula (14. 1 1 1) where 5 = ^din case (i) and 5 = din case (ii). 
By Remarks 12. 101 and 14.61 one concludes that the spectrum of the block operator matrix L con- 
sists of a two disjoint components Oq = spec(Zo) and g[ = spec(Zi) lying in the closed r- 
neighborhoods Or(spec(Ao)) and Or(spec(Ai)) of the corresponding spectral sets spec(Ao) and 
spec(Ai). 

Remark 4.14. Examples |4.15| and |4.16| below show that the bound ( 14.201 ) is sharp in the fol- 
lowing sense. Given a number d >0 and values of the norms ||B|| and ||C|| satisfying (14.191 ) one 
can always present self-adjoint (and even rank one or two) entries Aq, Ai and bounded B and C 
such that in case (ii) the bound (14.201 ) turns into equality. Notice that Examples 14.151 and 14.161 
serve for the spectral dispositions (14.131 ) and (14.141) . respectively. 

Example 4.15. Let = ^i = C. In this case the entries Aq, Ai, B and C of (12.11 ) are simply 
the operators of multiplication by numbers. Set Aq = — |, Ai = ^, B = b, and C = —c where 
b,c, and d are positive numbers such that \/^ < d/2. The Riccati equation (12.11 ) turns into a 
numeric quadratic equation whose solutions K^^^ and K^^^ are given by 

K^'^ = , T^P)^ £ (4.21) 

The right-hand sides of the equalities in (14.211) also represent the norms of the corresponding 



d 

2W 



solutions K^^^ and K^'^\ Obviously, only the solution K^^'' satisfies the bound HA"!! < 
Also notice that the eigenvalues of the associated 2x2 matrix L (which is given by ( 12.71 )) read 
A_ = — \Jd'^ jA — bc and A+ = — A_. One observes, in particular, that A_ = Aq^BK^^\ 

Example 4.16. Let i^o = C and 9)\ = C^. Assume that 

Ao = 0, M = ^\ 2), B=(0 b), and C = f ^ 



where b^c, and d are positive numbers such that \/bc < d/2. In this case the Riccati equation 
(12.11) is easily solved explicitly. It has two solutions K^^^ = f ^^j^ 1 and K^^^ = f 1 with 
k^^^ = k^^^ = and 

^(1) - ^ ^(2) ^ 

A._j_ , , A._|_ 
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Clearly, ||S|| = b, \\C\\ = c, and only the solution K^^'' belongs to the ball ^dl{2\\B\\){^\^^o)- Its 
norm is given by the equality 

IICII 



\K 



(1)1 



| + Vf-||fi||||c|| 



5. /-SYMMETRIC PERTURBATIONS 

In this section we deal with perturbations of spectral subspaces of a self-adjoint operator 
under off-diagonal 7-self-adjoint perturbations. 

For notational setup we adopt the following hypothesis. 

Hypothesis 5.1. Assume thatAo andAi are self-adjoint operators on the Hilbert spaces f^nd 
9)\ with domains Dom(Ao) andT)om[Ai), respectively. Let B be a bounded operator from to 
9)q and C = —B*. Also assume that A and V are operators on Sj = S)o(B Sji given by ( 12. 9b and 
(12.10b . respectively, and L = A + V with Dom(L) = Dom(A). 

By J, 

(cf. (11.3b ) we denote a natural involution on the Hilbert space ^ associated with its orthogonal 
decomposition = ijo © i ■ Subsequently introducing the indefinite inner product 

[x,y] = {Jx,y), x,yGS:, (5.2) 

turns S) into a Krein space that we denote by ^. 

A (closed) subspace £ C is called uniformly positive if there is 7 > such that 

[x,x] > 7 \\x\\^ for any nonzero x E ^. (5.3) 

The subspace £ is called maximal uniformly positive if it is not a subset of any other uniformly 
positive subspace of ^. Uniformly negative and maximal uniformly negative subspaces of K are 
defined in a similar way. The only difference is in the replacement of (15.31 ) by the inequality 
[x,x\ < — 7 that should also hold for all x G x 7^ 0. For more definitions related to the 
Krein spaces we refer to |!9l, |[30l . 

Clearly, under Hypothesis 15. II both V and L are 7-self-adjoint operators on S), that is, the 
products JV and JL are self-adjoint with respect to the initial inner product (•,•). This means 
that V and L are self-adjoint on the Krein space ^. 

The statement below provides us with a sufficient condition for a self-adjoint block operator 
matrix L on ^ to have purely real spectrum and to be similar to a self-adjoint operator on Sj. 
Notice that for the particular case where the spectra of the entries Aq and Ai are subordinated, 
say supspec(Ao) < infspec(Ai), closely related results may be found in Uj Theorem 4.1] and 
||371 Theorem 3.2]. 

Theorem 5.2. Assume Hypothesis 15.71 Suppose that the Riccati equation 

KAq-AiK + KBK = -B* (5.4) 
has a weak (and hence strong) strictly contractive solution K, \\K\\ < 1. Then: 

(i) The operator matrix L has a purely real spectrum and it is similar to a self-adjoint 
operator on Sj. In particular, the following equality holds: 

L = TAT \ (5.5) 
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where T is a bounded and boundedly invertible operator on Sj given by 

I K* \ f I-K*K ^ "'^^ 
K I [ I-KK* 



(5.6) 



and A is a block diagonal self-adjoint operator on f), 

A = diag(Ao,Ai), Dom(A) = Dom(Ao) eDom(Ai), (5.7) 
whose entries 

Ao = (/ - K*Ky/^ (Ao + BK) {I - K*K)-'I\ 

Dom(Ao) = Ran(/-i^*^)V2|^^^^^^^, ^^-^^ 

and 

Ai = {I - KK*Yl'^{Ai- B*K*){I - KK*y^l^ , 

Dom(Ai)=Ran(/-^^*)V2|^^^^^^^, (5-9) 

are self-adjoint operators on the corresponding component Hilbert spaces i^o and 
(ii) The graph subspaces S^'q = '^(K) and Sj\ = ^[K*) are invariant under L and mutually 
orthogonal with respect to the indefinite inner product (15.21 ). Moreover, ^ = 
where the sign " [+] " stands for the orthogonal sum in the sense of the Krein space ^ 
The subspace S^'q is maximal uniformly positive while maximal uniformly negative. 
The restrictions of L onto the subspaces iip ^'^'^ ^'i ^''^ ^-unitary equivalent to the 
self-adjoint operators Aq and Ai, respectively. 

Proof. In the case under consideration the second Riccati equation (12.111 ) associated with the 
operator matrix L reads 

K'Ai -AqK' -K'B*K' = B. (5.10) 

Thus, it simply coincides with the corresponding adjoint (12.61 ) of the Riccati equation (15.41 ). By 
Remark 122] this means that the adjoint of K, K' = K*,is a. weak (and hence strong) solution to 
(15.101 ). Since \\K* \\ = \\K\\ < 1, the operators / - K*K and / - KK* are strictly positive, 

I-K*K>I-\\Kf >0 and I - KK* > I - \\Kf > 0, (5.11) 



and, hence, boundedly invertible. This also means that the operator T in (15.61 ) is well defined 
and bounded. In addition, by Remark l277] this implies that the operator W in (12.201 ) is boundedly 
invertible and, consequently, the same holds for T. 
Now notice that by Lemma 1231 we have 

Rm{K*K\uom{Ao)) C Dom(Ao), Ran(^^*|Do^(^,)) C Dom(Ai), 

and 

{I-K*K){Aq + BK)x={Aq + K*B*){I-K*K)x forall ;cGDom(Ao), (5.12) 

{I-KK*){Ax-B*K*)y = {Ai-KB){I-KK*)y for all j G Dom(Ai), (5.13) 

from which one easily infers that both Aq and Ai are self-adjoint operators. 

By using (15.81 ) and (15.91 ) one expresses the operators Zq =Aq+BK, Dom(Z()) = Dom(Ao), and 
Zi = Ai —B*K*, Dom(Zo) = Dom(A), in terms of Aq and Ai. Then combining the expressions 
obtained with equality (12.241 ) from Corollary 12.91 we obtain formula (15.51 ). The similarity (15.51) 
means, in particular, that spec(L) is a Borel subset of R. This completes the proof of part (i). 

The /-orthogonality of the subspaces and is obvious since for any jc, j G of the form 

x = xo®Kxq, xq G i3o, wAy = K*yi®yi, ji G ^i, (5.14) 



18 



S. ALBEVERIO, A. K. MOTOVILOV, AND A. A. SHKALIKOV 



we have [x,j] = {Jx,y) = {xQ,K*yi) — {Kx(),y\) = 0. Thus, the fist two assertions of part (ii) 
follow from Corollary [19] (ii). On the other hand, (15.141 ) yields \\x\\^ < (1 + ||/<r|P)||xo|p and 
IblP < (1 + ll-^lP)ll3'i IP' ^rid, hence, combined with (15. Ill ), it implies [x,x] > y\\x\\^ and [y,y] < 
-r\\yf where 7 = (1 - + > 0. This means that Sj^ and Sj[ are maximal 

uniformly positive and maximal uniformly negative subspaces, respectively. 

Now introduce the operators Tq = Wo{I - K*K)-^/^ and Ti = Wi {1-KK*)-^I^ where Wq and 
W\ are given in ( 12.281 ) assuming that K' = K*. Taking into account ( 15.81 ) and (15.91 ). the identities 
(12^271) of Corollary [22] (ii) then imply 

Tq-^L\s^,To = Ao and Tf ^Ll^-y/i = Ai. (5.15) 

Clearly, Ran To = RanTi = [TqXo, Toyo] = {xo,yo) for any xo,yo G Sjo, and [Tixi , Tiyi] = 
— {x\ ,yi ) for any x\,y\ G . This means that both Tq: ^jq—^ and Ti : $j\ ^ Sj\ are ^-unitary 
operators. Therefore, equalities (15.151 ) prove the remaining statement of part (ii). 

The proof is complete. □ 

Remark 5.3. By equalities (15.81 ) and (15.91 ) the self-adjoint operators Aq and Ai are similar to 
the operators 

Zo=Ao + BK, Dom(Zo) =Dom(Ao), and Zi =Ai-B*K*, Dom(Zi) =Dom(Ai) (5.16) 
respectively, and, thus, 

spec (Ao ) = spec (Zo) and spec (Ai ) = spec (Zi). (5.17) 

Notice that identities (15.121 ) and (15.131 ) imply that the operators Zq and Zi are self-adjoint on 
the corresponding Hilbert spaces and i^i equipped with the new inner products (/o,^o)f,o = 
{{I-K*K)fo,go)^^ and {fugi)^, = respectively. 

Remark 5.4. The requirement \\K\\ < 1 is sharp in the following sense: If there is no strictly 
contractive solution to the Riccati equation (15.41 ) then the operator matrix L may not be similar 
to a self-adjoint operator at all. This is clearly seen from the simple example below. 

Example 5.5. Let = ^i = C. Set Aq = — ^, Ai = ^, and B = b where b and d are positive 
numbers such that b > ^. If > d/2, the Riccati equation (15.41) has two solutions X^^^ = 

4b ~^ W ^'^^ ^^^^ ~ lh~ ^ ~ ^- Both X^^^ and X^^^ are not strictly contractive since 
II = ||X(^' II = 1. At the same time the spectrum of the matrix L consists of the two complex 

eigenvalues Ai = i\J b^ — ^ and A2 = — /y''^^^-^. If b = j, the equation (15.41 ) has the only 
solution X = I. In this case the spectrum of the matrix L is real (it consists of the only point 
zero) but one easily verifies by inspection that the only eigenvalue of L has a nontrivial Jordan 
chain and, thus, L is not diagonalizable. Therefore, in both cases b > d/2 and b = d/2 the 
matrix L cannot be made similar to a self-adjoint operator 

The next assertion represents a quite elementary corollary to Theorem l5.2[ 

Lemma 5.6. Let the assumptions of Theorem 15.21 hold. Assume, in addition, that the spectra 
o'q = spec(Zo) and o[ = spec(Zi) of the operators Zq andZi given by (15.161 ) are disjoint, that is, 
(Tq n a( = 0. Then <Jq and o[ are complementary spectral subsets of the block operator matrix 
L, spec(L) = o'q U o[, and the graphs i^Q = ^{K) and ^\ = ^{K*) are the spectral subspaces 
associated with the subsets <Jq and respectively. 
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Proof. By the assumption the spectra spec(Ao) = spec(Zo) = cJq and spec(Ai) = spec(Zi) = o[ 
(see Remark l531 ) of the self-adjoint operators Aq and Ai given by (15.81 ). (15.91) are disjoint. Hence, 
the spectral projections Ea(o'o) and Ea(cj() of the self-adjoint diagonal block operator matrix 
A = diag(Ao,Ai) associated with its spectral subsets Oq and a| read simply as 

EaK) = (o o) ^"'^ EA(a()=([; ^) 

By Theorem 15.21 (i) the operator L is similar to the operator A. This means that the similar- 
ity transforms ^i{Oq) = T^\{Oq)T^^ and Ei(a|) = rEA(cTj)r^^ of the spectral projections 
Ea(cJo) arid Ea(cj|) with T given by ( 15.61 ) represent the corresponding spectral projections of L. 
One verifies by inspections that = Q'S(k) and El(cJj') = Qc^^k') where Qc^{^k) and Q'^[k') 

are given by (12.231) assuming that K' =K*. That is, El^Gq) and Ei(a( ) are the (oblique) projec- 
tions onto the graph subspaces "^{K) and "^{K*), respectively, which completes the proof. □ 

Remark 5.7. The spectral projections EL(aQ) = Q^ik) and Ei(a() = are orthogonal 

projections with respect to the Krein inner product (15.2b . 

From now on we will assume that the spectra of the entries Aq and Ai are disjoint and, thus, 
the sets Gq = spec(Ao) and Gi = spec (A i) appear to be complementary disjoint spectral subsets 
of the total self-adjoint operator A. In such a case for any bounded perturbation V satisfying the 
bound ||y II <d/2,d = dist(ao, CJi), the spectrum of the perturbed operator L=A + V consists of 
two disjoint subsets a,, and g[, lying in the closed ||y|| -neighborhoods 0\\v\\ (co) and 0\\v\\ {<^\) 
of the spectral sets Gq = spec(A()) and Gi = spec(Ai), respectively. One can think of the sets Gq 
and g[ as the result of the perturbation of the corresponding spectral sets Gq and Gi . 

Provided that the perturbation V is /-symmetric and ||y|| < d/2, Theorem [53] below gives 
sufficient a priori conditions for the perturbed operator L = A + V to remain similar to a self- 
adjoint operator. Hence, this theorem also gives sufficient conditions for the perturbed spectral 
sets Gq and g[ to remain on the real axis. Furthermore, the theorem presents the main result of 
the section giving for such V ana priori norm bound on variation of the spectral subspaces of A 
associated with the disjoint spectral subsets Gq and Gi. 

Theorem 5.8. Assume Hypothesis \5.1\ and choose one of the following: 

(i) Assume ( 14.121 ) and set 5 = ^d; 

(ii) Assume (14.131 ) or (14.141) and set 5 = d. 

Also suppose that 

||V||<|. (5.18) 

Then the spectrum of the operator L is purely real and consists of two disjoint components Gq 
and g[ such that 

00 C 0,.(spec(Ao)) and a( C 6>r(spec(Ai)) , (5.19) 

where 

(\ 2||y||\ 

r = ||V||tanh -arctanh — - — < ||V||. 
V2 b ) 

Moreover, the operator L is similar to a self-adjoint operator and the same is true for the parts 
ofL associated with the spectral subsets Gq and g[. Furthermore, the following bound holds: 

/\ 2||y||\ 

tanOo < tanh ( - arctanh — - — j , (5.20) 
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where ©o = 0(^O)^o) denotes the operator angle between the subspace and the spectral 
subspace 9)'q ofL associated with the spectral subset Oq. Exactly the same bound holds for the 
operator angle ©i = between the subspace and the spectral subspace of L 

associated with the spectral subset o[. 

Proof. Under either assumption (i) or (ii) from Theorem l4.11l it follows that the Riccati equation 
(15.41 ) associated with the block operator matrix L has a solution K G ^(^0,^1) that is unique in 
the ball ^s/2\\b\\ (^Oj^i) and satisfies the bound (see formulas (14.181 ) and (14.201 )) 



Here we have taken into account that ||B|| = ||y ||. We refer to Remark |43] regarding the use of 
the hyperbolic tangent in (15.211) . 

Clearly, the bound (15.211) yields that the solution Kisa strict contraction, < 1. Then by 
Theorem l5.2l the block operator matrix L is similar to the self-adjoint operator A given by (15.71) - 
(15.91) . Hence spec(L) C M and spec(L) = o'^VJ o[ where a,, = spec(Ao) and o[ = spec(Ai ). By 
Remark [531 we also have a,, = spec(Zo) and o[ = spec(Zi) where Zo and Z\ are given by (15.161 ). 
Since \\BK\\ < ||V||||/r|| < r and \\B*K*\\ < ||V||||^|| < r, for the spectral sets o'q = spec(Zo) 
and o[ = spec(Zi) the inclusions (15.191 ) hold and these sets are disjoint, dist(aQ, a() > 5 — 2r> 
5 — 2||y|| >0. To prove the remaining statements of the theorem one only needs to apply Lemma 
15.61 and then to notice that due to (12.141) we have || tan ©oil = ||^|| and hence tan©o < \\K\\. 
Similarly, tan©i < \\K*\\ = \\K\\. 

The proof is complete. □ 

Remark 5.9. By the upper continuity of the spectrum, the inclusion spec(L) C M also holds for 
||V|| = d/nm case (i) and for ||y || = d/2 in case (ii). 

Remark 5.10. In case (ii) the bounds ( 15. 19b on the location of spec(L) and the bound ( |5.20b on 
the angle ©0 are optimal. The optimality of both ( 15. 19b and ( 15.201 ) is seen from Examples 14.151 
and |4. 161 where one sets c = b. 

Remark 5.11. Under condition (15.18b in both cases (i) and (ii) the perturbed spectral subspaces 
and are mutually orthogonal with respect to the Krein space inner product (15.2b and, 
thus, R = ^q[+]^j. These subspaces are maximal uniformly positive and maximal uniformly 
negative, respectively. The restrictions of L onto S)'q and ^\ are .^-unitary equivalent to the self- 
adjoint operators Aq and Ai given by (15.81 ) and ( 15.91 ). respectively. By Theorem 15.21 (ii) all this 
follows from the fact that \\K\\ < 1 which we established in the proof of Theorem 15. 81 

Theorem 15 . 8 I claims that the spectrum of the block operator matrix L is purely real whenever 
the off-diagonal 7-self-adjoint perturbation V satisfies the bounds ||V|| < d/2 in case (i) or 
\\V\\ < d/n'm case (ii). Recall that case (ii) corresponds to the general spectral situation where 
no constraints are imposed on the mutual positions of the spectra spec(Ao) and spec (A 1) except 
for the condition (14.121 ). Now we want to prove that, in fact, under the only condition (14.12b the 
spectrum of the operator L remains purely real even if d/n < \\V\\ < d/2, at least in the case 
where the entries Aq and Ai are bounded. Our proof will be based on results from ll32l and H9l . 

Theorem 5.12. Assume Hypothesis \5.1\ Assume, in addition, that both the entries Aq and Ai 
are bounded and such that dist(spec(Ao),spec(Ai)) = d > 0. Also suppose that \\V\\ < d/2. 
Then the spectrum of the block operator matrix L is real, that is, spec(L) C M. 




(5.21) 
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Proof. Under Hypothesis 14. 8 1 and condition \\V\\ <d/2 the inclusion spec(L) C M has been al- 
ready proven in Theorem 15.81 (ii). Thus, let us only consider the case that is not covered by 
Hypothesis 14.81 In this case, because of the separation condition dist(spec(A()),spec(Ai)) = d, 
the specti^um of Aq consists of several (at least two) nonempty subsets isolated from each 
other at least by the distance Id. Denote these isolated spectral subsets of Aq by Oq \ i = 
1,2,. . . ,no,no>2, assuming that they are numbered from left to right (i.e. supa^'^ < inf a^'^''), 
the gap between supa^'^ and info^'^^^ contains a nonempty subset of the spectrum of Ai, 
and |J"=i f^d'^ = spec(Ao). In exactly the same way, divide the spectrum of Ai into the sub- 
sets aP, j = l,2,...,ni, ni > 2, so that \J%\(^i^^ = spec(Ai), supa^ < infaj'-'^^^ and 

(sup(T|^\infaj^^+'^)nspec(Ao) / 0. Denote by ^d'^, /= 1,2,..., no, and Sj[^\ j = l,2,...,nu 
the spectral subspaces of the operators Aq and Ai associated with the corresponding spectral 

subsets ad'' and o[^K Surely, eZa^o^ = and = ^H- 

Now take arbitrary unit vectors 

4''^e-^o\ ||4'''|| = 1, /=l,2,...,?io, and G \\e['^ \\ = 1, j = 1,2, . . . ,ni, (5.22) 
and construct numerical matrices Aq, Ai, and B with the entries 

Aojt = {Aoef\ef), A^j, = {A,e\'\ef), and B^j = {Bef\e'^''^), 

respectively. Consider the matrices Aq and Ai as operators resp. on = C"" and = C"', 
and B as an operator from S)i to i^o- Out of the matrices Aq and Ai construct the block diagonal 

matrix A = diag(Ao, Ai) and out of B and B* the off-diagonal matrix \/ = ^ ) Both 

matrices A and V have dimension n xn where n = no + ni, and we consider them as operators 
on the n-dimensional space S) = S^o® ^i- 

Our nearest goal is to prove that the spectrum of the operator L = A + V is real. To this end, 
first, introduce the indefinite inner product 

[x,y] = {xQ,yo)^^^-{xi,yi)^^, x = xo®xi,y = yoQyi, xo,yo G^c-^bJi G^i, (5-23) 

which turns the Hilbert space Sj into a Krein (Pontrjagin) space. We denote the latter by K The 
operator A is self-adjoint both on and ^ while B only on ^. 

Then notice that for different / and k the vectors ef ' and ef'^ belong to the different (and mutu- 
ally orthogonal) spectral subspaces of Aq and, hence, Ao,;^ = ^f^5ik where Xf^ = {Aoef'\ef^) 
and Sjji is the Kronecker's delta. Similarly, Aij-/ = Aj^'Sj; where Aj^' = {A\eJ\e'j'^). Clearly, 
both Xf\ i = 1,2,. . . ,no, and j = 1,2,... ,n\, are simple eigenvalues of A and, by con- 
struction of A, one has Xf^^ G cony{alf') and Aj'"* G conv(aj^'). This yields 

min |A/°'-Af'| >2J, min lAp' - > 2J, and minlA^' - > t/. (5.24) 

It is also obvious that, with respect to the inner product (15.231 ). the eigenvalues A/°' , / = 1 , 2, . . . , ?io, 
are of positive type, while the eigenvalues Xj^\j = 1,2, . . . are of negative type. 

Now to prove the inclusion spec(L) C M it only remains to observe that ||V|| < ||y|| < d/2 
and then to apply ||32l Corollary 3.4] (cf. HH Theorem 1.2]). 
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Since the inclusion spec(L) C M holds for any choice of the vectors (15.22b . one then concludes 
that also W"{L) C M where VK"(L) denotes the block numerical range (see B9l Definition 2.1]) 
of the operator L with respect to the decomposition 



By ||49l Theorem 2.5] we have spec(L) c W"{L). Hence, spec(L) c M, which completes the 
proof. □ 

Remark 5.13. By the upper continuity of the spectrum, under the hypothesis of Theorem 15.121 
the spectrum of L = A + V is real also for ||y || =d/2 (cf. Remark ls!9l ). 

Remark 5.14. Under the assumptions of Theorems 15 . 8 1 (ii) or l5.12l the requirement ||y|| <d/2 
guaranteeing the inclusion spec(L) C M is sharp. This is seen from Example [53] with b > d/2. 

6. Quantum harmonic oscillator under a ^^^-symmetric perturbation 

Let A be the Schrodinger operator for a one-dimensional quantum harmonic oscillator (see, 
e.g., 138] Chapter 12]). The corresponding Hilbert space is = LjiM.). Assuming that the units 
are chosen in such a way that h = m = co = I, the operator A reads 

(A/)W = -^^/(x) + i^VW, Dom{A) = !^fGWi{R) J^dx ix)\^ Koo^ (6.1) 

where W2^(M) denotes the Sobolev space of those L2(M) -functions that have their second deriva- 
tives in L2(M). The subspaces 

^0 = Z.2,even(ffi) and fli = L2,odd(M) (6.2) 

of even and odd functions are the spectral subspaces of the (self-adjoint) operator A associated 
with the spectral subsets 

Go = spec (A 1^^^) = {n + \/2 | n = 0,2,4, ... } and ai = spec (A j^^^) = {?i+l/2 | n= 1,3,5...}, 

respectively (see, e.g., Ii44l p. 142]). Clearly, 9j = ^o©^i> the spectral sets Oq and Oi are 
disjoint, 

d = dist{ao,ai) = I, and aoU ai = spec(A). (6.3) 

Let be the parity operator on L2(M), {j^f){—x) = f{—x), and J7 the (antilinear) operator 
of complex conjugation, {£^f){x) = f{x), f G L2(M). An operator V on L2(M) is called 3^ S^- 
symmetric if it commutes with the product ^ S/' , that is, 3"V = V.^^ (see, e.g. llT5l[T6l 
and references therein). 

In a particular case where the ^^^-symmetric potential V is an operator of multiplication by 
a function V{-) of Loo(M), the following equality holds (see, e.g., I3j|; cf. [i32il ): 



V{x) = V{-x) for a.e. xeR (6.4) 

and hence 

V* = m^. (6.5) 

Observe that the parity operator represents nothing but the involution ( 11.31 ) associated with 
the complementary spectral subspaces (16.2b of the oscillator Hamiltonian (16.11 ). Therefore, the 
equality (16.51 ) implies that the 3^ ^ -%^mmsXnc multiplication operator W is 7-self-adjoint with 
respect the involution J = ^ . 
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Any bounded complex-valued function V on M possessing the property (16.41) admits the rep- 
resentation 

V{x)=a{x) + ib{x) (6.6) 
where both a and b are real-valued functions such that 

a{—x) = a(x) and b{—x) = —b{x) for any a; G M. 

The terms Vdiag (x) = a{x) and Voff{x) = ib{x) represent the corresponding parts of the multiplica- 
tion operator V that are diagonal and off-diagonal with respect to the orthogonal decomposition 
^ = -Qo ©^1, that is, with respect to the decomposition L2(M) = L2,even(IR) ©^2,odd(IR)- 

Now assume that V is an arbitrary bounded off-diagonal operator on ^ = L2 (M) being 7-self- 
adjoint with respect to the involution J = 3^ . One can choose in particular a ^j7-symmetric 
potential (16.61 ) with a = 0. By taking into account (I6.3I ). from [ 15 , Theorem 1.2] it follows that 
the spectrum of the perturbed oscillator Hamiltonian L = A + V, Dom(L) = Dom(A), remains 
real (and discrete) whenever ||y|| < 1/2. If, in addition, the bound ||y|| < I/tt is satisfied then 
one can tell much more: Under such a bound Theorem 15.81 (i) implies that L is similar to a 
self-adjoint operator. This theorem also gives bounds on the variation of the spectral subspaces 
(lO) : 

tanQy < tanh Qarctanh(7r||y||)^ < 1, j = 0,l, 

where = stands for the operator angle between the subspace and the spectral 

subspace of the perturbed oscillator Hamiltonian L = A + V associated with the spectral 
subset o'j = spec(L) n 0\\v\\ (ay), j = 0, 1. 

Acknowledgments. The authors thank S.M. Fei for his useful remarks on ^j7-symmetric 
operators. A. K. Motovilov and A. A. Shkalikov gratefully acknowledge the kind hospitality of 
the Institut fiir Angewandte Mathematik, Universitat Bonn, where the main part of this research 
has been performed. 



References 

[1] V. M. Adamjan and H. Langer, Spectral properties of rational operator valued functions, J. Open Th. 33 (1995), 
259 - 277. 

[2] V. Adamyan, H. Langer, and C. Tretter, Existence and uniqueness of contractive solutions of some Riccati 

equations, J. Funct. Anal. 179 (2001), 448 - 473. 
[3] S. Albeverio, S. M. Fei, and P. Kurasov, Point intereactions: £P £^-Hermiticity and reality of the spectrum, Lett. 

Math. Phys. 59 (2002), 227 - 242; arXiv: quant-ph/02061 12 
[4] S. Albeverio and S. Kuzhel, Pseudo-hermicity and theory of singular perturbations, Lett. Math. Phys. 67 

(2004), 223 - 238. 

[5] S. Albeverio, K. A. Makarov, and A. K. Motovilov, Graph subspaces and the spectral shift function, Canad. J. 

Math., 55 (2003), 449 - 503; math.SP/0105142 v3. 
[6] S. Albeverio and A. K. Motovilov, Operator integrals with respect to a spectral measure and solutions to some 

operator equations. Fundamental and Applied Mathematics (to appear); arXiv: math.SP/0410577 v2. 
[7] S. Albeverio, A. K. Motovilov, and A.V. Selin, The a priori tanS theorem for eigenvectors, SIAM J. Matrix 

Anal. Appl. 29 (2007), 685 - 697; math.SP/05 12545 
[8] W. Arendt, F. Rabiger, and A. Sourour, Spectral properties of the operator equation AX + BX — Y, Quart. J. 

Math. Oxford 45 (1994), 133 ~ 149. 
[9] T. Y. Azizov and I. S. lokhvidov. Linear Operators in Spaces with an Indefinite Metric, John Wiley & Sons, 

Chichester, 1989. 

[10] C. M. Bender, Making sense of non-Hermitian Hamiltonians, Rep. Prog. Phys. 70 (2007), 947 - 1018; arXiv: 
[hep^th/0703096| 



24 



S. ALBEVERIO, A. K. MOTOVILOV, AND A. A. SHKALIKOV 



[11] C. M. Bender and S. Boettcher, Real spectra in non-Hermitian Hamiltonians having IPS/' Symmetry, Phys. 

Rev. Lett. 80 (1998), 5243 - 5246; arXiv:|i physics/9712001| 
[12] C. M. Bender, S. Boettcher, and P. N. Meisinger, S/" -symmetric quantum mechanics, J. Math. Phys. 40 

(1999), 2201 - 2229; arXiv: quant-ph/9809072 
[13] R. Bhatia, C. Davis, and A. Mcintosh, Perturbation of spectral subspaces and solution of linear operator 

equations. Linear Algebra Appl. 52/53 (1983), 45 - 67. 
[14] R. Bhatia and P. Rosenthal, How and why to solve the operator equation AX — XB = Y , Bull. London Math. 

Soc. 29(1997), 1 -21. 

[15] E. Caliceti, F. Cannata, and S. Graffi, Perturbation theory ofPT symmetric Hamiltonians, J. Phys. A 39 (2006), 

10019 - 10027; arXiv: math-ph/0607039 
[16] E. Caliceti, S. Graffi, and J. Sjostrand, Spectra of PT -symmetric operators and perturbation theory, J. Phys. A 

38 (2005), 185 - 193; arXiv: math-ph/0407052 
[17] C. Davis and W. M. Kahan, The rotation of eigenvectors by a perturbation. Ill, SIAM J. Numer. Anal. 7 (1970), 

1 -46. 

[18] J. W. Demmel, Three methods for refining estimates of invariant subspaces. Computing 38 (1987), 43 - 57. 

[19] K.-J. Engel and R. Nagel, One-Parameter Semigroups for Linear Evolution equations. Graduate Texts in Math- 
ematics 194, Springer- Verlag, New York, 2000. 

[20] M. Haase, The Functional Calculus for Sectorial Operators and Similarity Methods, Dr. rer. nat. thesis, Uni- 
versitiit Ulm, 2003. 

[21] P.R. Halmos, Two subspaces. Trans. Amer. Math. Soc. 144 (1969), 381 - 389. 

[22] V. Hardt, A. Konstantinov, and R. Mennicken, On the spectrum of product of closed operators. Math. Nachr. 
215 (2000), 91 - 102. 

[23] E. Heinz, Beitrdge zur Storungstheorie der Spektralzerlegung, Math. Annalen 123 (1951), 415 - 438. 
[24] T. Kato, Perturbation Theory for Linear Operators, Springer- Verlag, Berlin, 1966. 

[25] V. Kostrykin, K. A. Makarov, and A. K. Motovilov, Existence and uniqueness of solutions to the operator 
Riccati equation. A geometric approach. Contemporary Mathematics (AMS) 327 (2003), 181 - 198; arXiv: 
math.SP/0207125 

[26] V. Kostrykin, K. A. Makarov, and A. K. Motovilov, On the existence of solutions to the operator Riccati equa- 
tion and the tan® theorem, Integr. Eq. Oper. Th. 51 (2005), 121 - 140; arXiv: math.SP/0210032 v2. 

[27] V. Kostrykin, K. A. Makarov, and A. K. Motovilov, Perturbation of spectra and spectral subspaces. Trans. 
Amer. Math. Soc. 359 (2007), 77 - 89; arXiv: math.SP/0306025 

[28] V. Kostrykin, K. A. Makarov, and A. K. Motovilov, A generalization of the tan 20 Theorem, Operator Theory: 
Adv. Appl. 149 (2004), 349 - 372; arXiv: math.SP/0302020 

[29] D. Krejcifik, Calculation of the metric in the Hilbert space of a ^ .S/' -symmetric model via the spectral theorem, 
J. Phys. A 41 (2008), 244012 (6 pp.); arXiv:0707.1781 

[30] H. Langer, Krein space, in: Encyclopaedia of Mathematics (Ed. M. Hazewinkel); |http://eom.springer.de^ 
k055840.htm. 

[31] H. Langer, A. Markus, V. Matsaev, and C. Tretter, A new concept for block operator matrices: the quadratic 
numerical range. Linear Algebra Appl. 330 (2001), 89-1 12. 

[32] H. Langer and C. Tretter, A Krein space approach to PT-symmetry, Czech. J. Phys. 54 (2004), 1113 - 1120; 
Corrigendum, Ibid. 56 (2006), 1063 - 1064. 

[33] H. Langer and C. Tretter, Diagonalization of certain block operator matrices and applications to Dirac opera- 
tors. Operator Theory: Adv Appl. 122 (2001), pp. 331 - 358. 

[34] G. Lumer and M. Rosenblum, Linear operator equations, Proc. Amer. Math. Soc. 10 (1959), 32-41. 

[35] R. McEachin, Closing the gap in a subspace perturbation bound. Linear Algebra Appl. 180 (1993), 7 - 15. 

[36] R. Mennicken and A. K. Motovilov, Operator interpretation of resonances arising in spectral problems for 
2x2 operator matrices. Math. Nachr. 201 (1999), 1 17 - 181; arXiv: funct-an/970800_l 

[37] R. Mennicken and A. A. Shkalikov, Spectral decomposition of symmetric operator matrices. Math. Nachr. 179 
(1996), 259-273. 

[38] A. Messiah, Quantum Mechanics, Vol. /, Wiley & Sons, 1963. 

[39] A. Mostafazadeh, Pseudo-Hermiticity versus PT symmetry: The necessary condition for the reality of the 
spectrum of a non-Hermitian Hamiltonian, J. Math. Phys. 43 (2002), 205 - 214; arXiv: math-ph/0 107001 

[40] A. K. Motovilov, Removal of the resolvent-like energy dependence from interactions and invariant subspaces 
of a total Hamiltonian, J. Math. Phys. 36 (1995), 6647 - 6664; arXiv: ,funct-an/9606002i 



BOUNDS ON VARIATION OF SPECTRAL SUBSPACES UNDER /-SELF-ADJOINT PERTURBATIONS 



25 



[41] A. K. Motovilov and A. V. Selin, Some sharp norm estimates in the subspace perturbation problem, Integr. Eq. 

Oper. Th. 56 (2006), 5 1 1 - 542; arXiv: |math.SP/0409558| v2. 
[42] M. T. Nair, An iterative procedure for solving the Riccati equation AiR — RAi = A3 +RA4R, Studia Math. 147 

(2001), 15-26. 

[43] V. Q. Phong, The operator equation AX — XB = C with unbounded operators A and B and related abstract 

Cauchy problems. Math. Z. 208 (1991), 567 - 588. 
[44] M. Reed and B. Simon, Method of Modern Mathematical Physics, I: Functional Analysis, Academic Press, 

1980. 

[45] M. Rosenblum, On the operator equation BX - XA = Q, Duke Math. J. 23 (1956), 263 - 269. 

[46] G. W. Stewart, Error and perturbation bounds for subspaces associated with certain eigenvalue problems, 

SIAM Review 15 (1973), 727 - 764. 
[47] G. W. Stewart, Error bounds for approximate invariant subspaces of closed linear operators, SIAM J. Numer. 

Anal. 8(1971), 796-808. 

[48] T. Tanaka, General aspects of PT -symmetric and P-self-adjoint quantum theory in a Krein space, J. Phys. A 39 

(2006), 14175 - 14203; arXiv: hep-th/0605035 
[49] C. Tretter and M. Wagenhoffer, The block numerical range of an nx n block operator matrix, SIAM J. Matrix 

Anal. Appl. 24 (2003), 1003 - 1017. 
[50] M. Znojil, Solvable PT-symmetric Hamiltonians , Phys. Atom. Nucl. 65 (2002), 1149 - 1151; 

arXiv: quant- ph/0008125 

Sergio Albeverio, Institut fiir Angewandte Mathematik, Universitat Bonn, Wegelerstrasse 
6, D-53115 Bonn, Germany; SFB 611 and HCM, Bonn; BiBoS, Bielefeld-Bonn; CERFIM, Locar- 
no; Accademia di Architettura, USI, Mendrisio 

E-mail address: albeverioSuni-bonn . de 

Alexander K. Motovilov, Bogoliubov Laboratory of Theoretical Physics, JINR, Joliot-Cu- 
RiE 6, 141980 Dubna, Moscow Region, Russia 
E-mail address: motovilvOtheor . j inr . ru 

Andrei A. Shkalikov, Faculty of Mathematics and Mechanics, Moscow Lomonosov State 
University, Leninskie Gory, 1 19992 Moscow, Russia 
E-mail address: ashkalikovSyahoo . com 



